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Abstract 

Using a two-dipole model of an optical near-field of Scanning Near-field Optical Microscope tip, 
i. e. taking into account contributions of magnetic and electric dipoles, we propose and analyze a 
new type of 3D optical nanotrap found for certain relations between electric and magnetic dipoles. 
Electric field attains a minimum value in vacuum in the vicinity of the tip and hence such a trap 
is quite suitable for manipulations with cold atoms. 



*Electronic address: vklim@sci.lebedev.ru 



1 



Recent enormous progress in the study of laser-cooled atoms and molecules put forward 
the problem of their using for quantum computing, frequency standard construction, and 
other technological applications (see e.g. [1-3]). To achieve this goal, as well as for the further 
progress of fundamental experiments in the field, new compact traps and new methods of 
"handling" the cold atoms, for example their transportation to/from a trap or between 
different traps, should be elaborated. An example of successful work in this direction is 
given e.g. by "atom-chip technology" experiments [4, 5], where guiding of cold atoms along 
the wires has been demonstrated. Further miniaturization of such devices is desirable. 
Ideally, one should have at his/her disposal true "single cold atom nanotweezers" which 
have very small sizes and only slightly perturb the trap. Evidently, the near-field optical 
configurations, based on subwavelength-size aperture in the apex of a sharp fiber tip or local 
field enhancement in the vicinity of sharp conducting tips (see e.g. [6, 7] for review on near- 
field optics) look rather promising. Hence it is not surprising that a number of near-field 
optical traps/tweezers have been proposed [7 - 15]. 

However, all of them have certain drawbacks and, we believe, this is the main reason 
why, to the best of our knowledge, none has been realized up to now. Putting aside the 
configurations where an extremum of the optical field is achieved on the surface of the tip 
[13] (configuration obviously inappropriate for atom trap), among these drawbacks we could 
mention that the traps proposed needed different hard-to-control non-optical interactions 
(centripetal potential, gravitational interaction, van der Waals force, etc.) to be closed, can 
be realized only in the non practical light reflection mode from the subwavelength aperture, 
and so on. 

Here we propose and analyze the true "free standing" , or "support-free" purely optical 3D 
trap emerging in vacuum in the vicinity of an aperture of the Scanning Near-field Optical 
Microscope (SNOM) tip. The characteristic size of our trap is small in comparison with 
laser wavelength , so one indeed can speak about nanotrap. Because methods of Angstrom- 
precision motion of SNOM tip are well elaborated, the same construction is cold atom 
nanotweezers. 

Our analysis is based on the two-dipole model of optical near field occurring in the vicinity 
of this tip, see Fig. ^ Such a model has been established recently, when it has been shown 
that in addition to the classic Bethe and Bouwkamp consideration, where in the case of 
normal incidence optical near field is modeled by one magnetic dipole M (see e.g. Refs. 
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[16-18]), a field of an electric dipole P should be added to describe correctly optical near- 
field of a real fiber tip ^|-|^|. This effect is due mainly to the conical shape of the end 
area of such tip, and for a variety of tips and light polarizations different relations between 
the values and mutual orientations of M, P dipoles can be anticipated. 

In near field of tip as in any laser fields, resonant atoms are subject to the optical dipole 



force with the potential 



22 



^T ln(1 + W )f *lT^ (1) 

Here fi is an atomic transition dipole moment, 7 2 = Q 2 + (r/2) 2 where T is a natural 
line width, and Q = u — ujq is detuning between laser frequency u and resonant frequency 
of an atom ujq. Throughout the paper we consider the blue detuning fl > 0, which results 
in atom trapping at the minimum of laser field intensity.lt is well known that blue-detuned 
gradient force optical traps based on minimum of an electric field have essential advantages 
(they do not heat trapped atoms, etc.) in comparison with those based on maximum of an 



electric field using red-detuned light |23 1 . 



To characterize the proposed 3D nanotrap let us consider electric field of light in the 
vicinity of a SNOM tip (i. e. in the near-field regions of both dipoles). Corresponding 
electric field of an electric dipole P has the form: 

P R(RP) 

where R is the radius vector from the dipole position to an observation point. (SGCE units 
are used throughout the paper). Vector potential of a magnetic dipole M has the following 
form 

A=Ega ( 3) 

and according to Faraday's law electric field of this magnetic dipole in the near-field region 
has the form 

E m = ik 1 -^ (4) 

lu 

where k = — is wavevector in free space. When two dipoles are located at the same point 

c 

{0,0,0}, intensity of total electric field E = E m + E e behind an aperture can be presented 



in the form 

E 2 = |E e + E m | 2 = 

(5) 

k 2 \M\ 2 R 2 - k 2 (MR) (M*R) + |P| 2 + 3 (PR) (P*R) / R 2 - ik~R ([P*M] - [PM*]) 

R? 

Below we will speak about this value as about an "intensity" having in mind its obvious 
connection with the intensity of laser light Ii as "seeping" through the aperture of a SNOM 
tip: I^s — —E 2 , where c is speed of light . 

It is naturally to look first for an axially symmetric trap. Such a trap occurs in particularly 

if 

M z = P z = 
M x M* y + M y M* = 

\M X \ = \M y \ = M (6) 

p x p; + p v p* = o 

\Px\ = \Py\ = Pq 

If these conditions hold true, expression for the intensity can be rewritten in the form 

p2 2k 2 M 2 R 2 - k 2 M 2 p 2 + 2P 2 + 3P 2 p 2 /R 2 - zkQ 

where Q = i ([P*M] — [PM*]J and p = \J 'x 2 + y 2 is the distance from the symmetry axis 
to an observation point. 

Analysis shows that if we put the following additional condition on the dipole momenta 

Q 2 = 16M 2 P 2 (8) 
intensity (0) becomes equal to zero at the points 

x = y = 0, kz = — -r- I or [ x = y = 0, kz = — — ^- ) (9) 
In its turn, it is possible to show that (JSj) will be satisfied in the case 



P = P {i,l,0}, 
M = ±M o {z,l,0} 



(10) 
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It means that dipole momenta should be collinear to ensure the 3D trap. 

As intensity is a positive function of coordinates, zeros of intensity Q correspond to true 
3D minimum of an electric field. In FigJ2]the distribution of intensity in x — z plane is shown 
for the case P = 20; M = 10. 

Intensity J7J) has another point of extremum 

3 P \ ( 3 P \ 

x = 0, y = 0, kz = -—— ) or ( x = 0, y = 0, kz = — — I (11) 

y ' 2M J V 2M / v ; 

which is a saddle point (see Fig. 2.). The value of intensity at this point is 

E 2 = —k 6 ^l (12) 
729 P 4 V ; 

This quantity can be used to estimate the potential well depth. For such an estimation 
one can use Mq = Pq = Eqci 3 , where Eqis an amplitude of the incoming light wave in aperture 
plane and a is the radius of an aperture [16-19]. Substituting these values into ((HJ), we get 
for the potential well depth 

AE 2 ~— (kafEl (13) 
729 v ) o \ ) 

For modern tips ka ~ 1/2 and hence AE 2 ~ lO -3 ^. For typical experimental conditions, 
an intensity of optical near field at the aperture of SNOM tip is about I = 10 3 - lQ^W/cm 2 
[6,7]. This means that e. g. for alkaline atoms, which are characterized by resonant dipole 
moments of the order of 10~ 17 SGSE and natural line widths T = 10 7 — 10 8 s _1 (for example, 
for the 6S1/2-6P3/2 D2 transition of cesium atom at A=852 nm, /i=8.01-10 -18 CGSE and 
T = 3.07 ■ 10 7 s _1 [24]), traps with the depth of the order of a few milliKelvin, what is 
quite standard for optical dipole gradient force - based traps, can be realized when using 
Q = 100r, see (0). It is very important that under these conditions the trap will have about 
10 energy levels of atomic motion with lowest level being about lO -4 ^. 

It is worthwhile to note, that existence of such a 3D trap is highly nontrivial, because the 
relevant field components have very complicated structure (see Figs. EH 1110), an d it seems 
very difficult to provide a minimum for their sum. 

Mathematically, by varying P and M ratio and other parameters of the problem, our trap 
can be placed at any point on the symmetry axis. However, the two dipole approximation 
used is not valid very close to the aperture plane z = 0. Its validity starts from z ~ a, 
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where a is the radius of an aperture. Hence we should consider only such parameters where 
the trap position occurs not too close to the aperture plane. Besides the condition | z |> a 
allows us to neglect van der Waals attractive force which is always important in close vicinity 
of tip surface. 

For our trap the minimum of intensity is stable against small perturbations. For example, 
if relation between the dipole values is 

P= ftp, 1,0}, 

M = ±M e i f{i, 1,0} 

then true 3D minimum still exists provided sin |<^| < | or — 0.2 < 99 < 0.2. Small variations 
of z-components of momenta (P z and M z ) also result in small variations of trap position. 
Small mutual displacements (splitting) of dipoles in axial and/or radial directions also result 
in variation of the trap. In the case when magnetic dipole is moving away from the trap, 
that is in the case when electric dipole is placed between the trap and the magnetic one, 
the trap suffers only minor shifts. In the case when magnetic dipole is placed between an 
electric one and the trap, the intensity minimum disappears for large enough splitting, see 
Fig© 

The retardation effects in the near field region generally are small. Nevertheless to esti- 
mate their influence let us consider the full electric field of electric and magnetic dipoles: 



r . 1 ik k 2 \ R(RP) / 3 3ik k 2 
E e = P — =r + -ttt + 1 ' 



R? R? R ) R? \R? R 2 R 
[MR] ( ik k 2 

m ~ R \R 2 R 



(15) 



Analysis shows that for such a case 3d trap survives only when £ = P /M < 1.5035. 
Positions of the trapping region and the saddle point are given respectively by the formulae 



( x = y = 0, kz = £ - e + 3£ 5 + ...) or (x = y = 0, kz = - (f - £ 3 + 3£ 5 + ...)) (16) 

and 



x 



0, kz 



3 , 9 ^3 
2 s 16 



189 
128^ 



£ 5 + ... I or ( x = y = 0, kz 



_ ( V^3_18?, 
1 2 s 16 128 



(17) 
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Intensity at the bottom of the potential well is now nonzero 

E 2 min = 2k e MZ(e-e+-) (is) 

while the intensity at the saddle point is 

^e = ^M 2 (^ + ^ + ...) (19) 

These values determine the depth of our trap. Comparing results (|16 p -([19| h where re- 
tardation effects are taken into account, with the quasistatic results we see that in the 
case of magnetic dipole domination (£ < 1) the position of the trap remains in the near field 
region and hence the retardation effects have only minor influence. On the other hand, in 
the case of substantially large amplitude of electric dipole, the retardation effects destroy 
our trap. 

Hence we have shown that minimum of an electric field with the size smaller than the light 
wavelength and depth of a few milliKelvin did occur in vacuum in the vicinity of a SNOM 
tip. This minimum is stable against perturbations and can be tuned both in position and 
depth by changing relation between Pq and Mq and incident light electric field amplitude 
EO. This attests the trap proposed as a very promising basic element for future cold atom 
nanotraps and nanotweezers. Finally, we would like to note the following. Despite the two 
dipole model of optical near-field, whose using is inherent to obtain the reported results, 
nowadays seems is well established and supported by experiments [19-21], complete and 
rigorous analysis of all possible dipole configurations occurring for near field of a real tip is 
still lacking. In particular, the conditions to be imposed on experimental setup to obtain 
the trapping configuration of magnetic and electric dipoles also remain to be understood. 
Nevertheless, we believe that broad possibilities to vary SNOM tip shapes and coatings [6, 
7] as well as to vary other parameters (e. g. incoming light polarization) give enough hope 
for practical realization of such a configuration. Indeed, it can quite happen that this is 
already done at least for one of plenty of the tips demonstrated up to date. 
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FIG. 1: Problem geometry 
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FIG. 2: Isolines of the intensity eq. ® in the plane y = for Pq = 10; Mo = 10 (natural 
logarithmic scale). 



11 




-4 -2 2 4 

kx 



2 

FIG. 3: Isolines of |E X | in the vicinity of the 3d trap in the plane y = for Pq = 10; Mo = 10 
(natural logarithmic scale). 
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FIG. 4: Isolines of |Ey| in the vicinity of the 3d trap in the plane y = for Pq = 10; Mo = 10 
(natural logarithmic scale). 
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FIG. 6: Variation of the trap position as a function of axial splitting between the dipoles. 
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